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Abstract 

The spin - spin mass splitting of light, heavy and mixed mesons are described 
within a good accuracy in the potential model with screened potential. We conclude 
that the long - distance part of the potential cannot be pure scalar and that a vector 
- scalar mixture is favoured. With the same parameters which gives correct average 
mass spectrum excellent spin - spin splittings of heavy quarkonia is obtained. The 
results are obtained by going beyond usually used perturbation method, namely 
using configuration interaction approach. 
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1 Introduction 



The problem of hyperfine splitting in mesons still attract wide interest. It is widely 
acceptable that quark potential model gives a rather good description of spin-average mass 
spectrum of hadrons , considered as composite system of quarks . However, the question 
of explaining the influence of spin, namely spin-orbit ("fine"), spin-spin ("hyperfine") 
interaction, is not solved yet. The problem of mass splitting is due to spin structure and 
it is closely connected with the Lorentz - structure of the quark potential. These effects 
are far from being solved yet. Quite recently QCD motivated potential was applied to 
heavy quarkonia. Unfortunately authors of [0] restrited themselves to comparing the 
results to a single experimental value J/'^ — rjc mass difference. No attempt was made 
to calculate other quark-antiquark pairs. Spectroscopy of heavy mesons is studied also 
in paper [§. The authors calculated only the spectra-mass of heavy mesons B, Bs, D, 
Ds with including spin-spin interaction. In 0] the first and second-order calculations 
of masses are in good agreement with the experimental data except for the higher spin 
states. Besides there are several works, in which the authors take limited number of pairs 
(like uu) of particles for which splitting are calculated. Their results are quite good, 
but in work 0] authors limited themselves by considering only tt, p, K, K* mesons. In 
analogy Buchmuller and Tye also considers only two systems, namely cc and hb though 
in states 15* and 25*. They had obtaine good results for these particular pairs. There is 
another paper P] , in which both splitting and decay properties are studied. In this paper 
also good results for hyperfine splitting are obtained. But the authors calculate only few 
especial mesons {J/"^ — rjc). Also the authors at calculations guess that in the hyperfine 
splitting gives contribution only the one - gluon term. From our paper it is visible, that 
confinement also plays an essential role at calculations the hyperfine effects in mesons. 

The main problem of our work is to clarify some aspects of hyperfine interaction in 
the framework of configuration interaction approach (CI approximation, or CIA 0]). 

Following many authors we assume admixture vector-scalar potential (soft model). 
We consider vector and scalar parts of static potential [0 



V(r) = Vs(r) + Vv(r) (1) 



where 



r OTT/X ^ ^ DTT/X ^ ^ 

and e is the mixing constant. 

The Hamiltonian can be written as 

H = Ho + Hss (3) 

where for screened potential: 

ifo = :^V^-^+./^(l-e--) (4) 
2m r Dvr/i ^ ^ 

where m- is reduced mass of qq- system and, h = c = 1 units are used. In (|^) only Hss 
we take into account, because we calculate hyperfine splitting in S*- waves. Then all terms 
in which contain the orbital quantum number / are absent. 
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In the frame of Breit-Fermi approach the spin-spin interaction term is 

Hss = sIsIaVv (5) 

where S 1,2 are the spin of the particles. 
Now, we consider Shroedinger equation 

{Ho + Hss) ^ (T^) = (r ) (6) 

Here we suggest to use CI approach, which was previously very successfully applied 
in atomic physics [0. The essence of this approximation is that wave function (r) (r) is 
expanded in set of eigenf unctions of the Hamiltonian Hq, is 

n 

With substituting (^) in to (|^) and using eigenvalue E'^ , we obtain homogeneous 
system of linear equations for a„ 

a„ (^^ - ^) = - ^ a„ {ipm \Hss\ Vn) (8) 

n 

which have to be truncated for reasonable large n. In is the eigenvalue of the nonper- 
turbative Hamiltonian: 

H^Lfra = E^ifm (9) 

It is evident, that the solution of ® exists only if the determinant, that contains of 
the coefficients are equal to zero. The diagonalization of this determinant gives the values 
of energies we are looking for. This is a good method to finding the eigenvalues. This 
procedure goes far outside of perturbative method. 



2 Hyperfine splitting 

In this work we obtained hyperfine splitting for heavy, light and mixed mesons using the 
configuration interaction approach. For example, we calculate the hyperfine splitting in 
and mesons with using oscillator potential. 

Let us write down for this case Fermi-Breit equation for the two-quark system in the 
form 

{-^A + Ar^-'^ + Hss^'^{j)=E^>{r) (10) 

where 

H^ = -^A + Ar^ (11) 
Hss = sls"2AVv (12) 
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and (|T2D is a additional for non-perturbatived Hamiltonian. 

In this case vector and scalar part of the potential is equal accordingly: 

Vv = -— + eAr^ and Vs = {I - e) Ar\ 
r 

From one hand the oscillator potential not so bad in describing quarkonia data, on the 
other hand allowing to obtain analitic basic solutions. Hss is the additional term, which 
have to be taken into account in CIA. Also one - gluon exchange term we put in to CIA. 
Below for obtaining the reliable results we shall use more realistic screened potential and 
make the numerical evaluation of the appropriate matrix elements. But for a while let us 
remain at using oscillator basis functions. Substituting the expanding the wave function 
into ( [To|) and using the condition of ortonormality of the basic functions we obtain the 
algebraic system of equations for the coefficients a„. 

Two remarks are to be done in connection with i\i2^ . First we shall consider the 
contribution to spin-spin term of the potential consisting of both one-gluon and many- 
gluon exchanges (^. Therefore 

Hss = ^ Slsl {Ae7Tas6 (7) + 6A) (13) 

S1S2 = — I for psevdoscalar mesons and S1S2 = | for vector mesons. In the latest 
calculations we shall take the interaction as a screened potential and we shall obtain the 
basic functions numerically and evaluate (|^) with these functions. The second remark 
concerns to the meaning of parameter e. According to general point of view |1| only 
vector contribution is present in (|^). In order to avoid introducing additional parameters 
we shall consider the mixing parameter e to be same for many-gluon terms approximately 
equal to e = 0.5. 

If we restrict ourselves by one term in expansion (|^) then (|^) reduce to perturbation 
method. We present the wave function in (|l^) as two basic function of (0), 



[r) = aiipi^r) +a2ip2[r) , (14) 
immediately leads to a much better approximation. Namely one obtains for the energies 

E,,2 = + + ± yiE^, - E^2 + Hu + H22f + m,2H2„ (15) 

where Hnm = (v^m \Hss\ ^n)- 

Let us write for completness the first two terms of oscillator wave-functions 



7r 



1/4 



<^is (r) 



1/4 



2 

3' 



1 - -(3r^ \ e'^Yl 



where (3 = y Arrig . In this case the corresponding matrix elements should be have the 
form 
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22 





Si S2 ^ 


AA + 




7ri/4 




37ri/2 




1 


AA + 




37rV4 


ml 


vri/2 



(16) 



rr _rr _ '^C^sP'^' , Si S2 i ^VGasP^/' \ 



Substituting (|TB|) into (|T3p and using standart values of E'^ for oscillator potential 
one can obtain the values for hyperfine splitting which are given in the Table 1. The 
parameters are chosen to be A = 0.014 GeV^, e = 1, rrig = 1.5 GeV, as = 0.32 since 
exactly these parameters give the best results for oscillator potential in meson masses. 

Tablel. 





1 


2 


3 


4 


5 


6 


7 


8 


9 


AMexp 


15 


37 


41 


44 


46 


47 


48 


49 


50 


51 


117 


25 




39 


43 


46 


48 


49 


51 


52 


53 


95 



Certainly using oscillator potential we obtain only didactic value. Therefore next we 
turned to use more realistic potential, namely, screened potential. The concrete form of 
such screened potential was previously suggested in 0. We choose screened potential, 
because it gives excellent description of mass-spectrum in nonrelativistic potential models 
for heavy mesons. As Gerasimov pointed out ||^ the QCD calculations on lattice indicates 
that the spin-spin forces are rather short range. Exactly the screened potential satisfies 
this condition. In this case even the basic solutions for unperturbed Hamiltonian can not 
be found in analitic form. We found these solutions numerically and evaluated the matrix 
elements numerically too. Final results for hyperfine splitting for screened potential are 

2 

given in Tables 1,2. The following parameters were used in screened potential |^ = 0.224 
GeV^, /i = 0.054 GeV. All parameters was took from |TD[. Experimental values were 
taken from 



Our calculations shown that the next terms of CIA method give contribution of order 
of 10% for heavy mesons and 35% for light mesons. Let us stress that the first term of 
CIA method is in fact just a perturbation result. CIA expansion better takes into account 
the interaction between particles. A similar approach was suggested in paper [l^, where 
expansion was carried in basic function of oscillator potential. The suggested method is of 
considerable interest, since the perturbation method is still used as the practical method 
present §, §, 
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3 Discussion 

For reasons we present our data in 2 tables. We compare our results with results obtained 
in the works [Q, ||13[, [Q. In the paper [0 good results are obtained, but the authors 
introduced additional parameters ro. In the paper [1^ hyperfine splitting is calculated in 
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the first account of perturbation tlieory. Tliere are sfiown tliat tlie first order perturbations 
theory gives a good description of experimental data. The best results in paper [0 



are 
In 



obtained. But only one - gluon exchange is taken in spin-spin forces in the paper |T4 
all these works there is one common fault as in them are restricted to viewing limited of 
number mesons. 

Table 2. 







AMtheor 
{MeV) 


@ 

AMtheor 
{MeV) 




AMtheor 
[MeV) 


Our results 

AMtheor 
[MeV) 


AMexp 
{MeV) 




634 


550 


651 


923 


635 




329 






411 


150 




217 




270 


580 


320 




135 






285 




AMk'-k 


405 


461 


393 


707 


398 




195 






336 


200 


AMd'-d 


92 


147 


150 


186 


143 


AMd"-d' 








112 




AMb'-b 


32 


52 


58 


57 


45.9 


AMb"-b' 








36 





Table 3. 





y 

AMtheor. 
{MeV) 


|13 

AMtheor. 
{MeV) 


y 

AMtheor 
{MeV) 


Our results 

AMtheor. 
{MeV) 


AMexp 
{MeV) 




87 


190 


128 


163 


144 










100 












50 


47 


AM„., „, 








33 












49 




AMb',_b'^ 








31 




AM^_^, 


31 


39 


82 


46 




AMy_,, 


9 






26 






65 


100 


112 


110 


117 


AM*./?;, 


32 


54 




67 


95 



In Table 2 we show final data for light-quark systems, which has mainly relativistic 
character and compare them with other data. The calculation are carried out by using 
(J^-^). In our quazirelativistic approach we should take into account in Hamiltonian also 
the term of order p^ i. In paper [14] it was shown that the spectrum of Hamiltonian 
in nonrelativistic potential models, and spectrum of relativistic Hamiltonian for bound 
state and for first radial excited state are equivalent. On the other hand Lucha and 
Shoeberl [|[] has shown that in some cases this relativistic kinematik term change the 
mass-spectrum drustically. But since we consider mass- difference therefore all relativistic 
effects must be cancel. In other papers authors introduce new additional parameters and 
obtain good description of hyperfine splitting. In paper Faustov et al. also had 
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made similar calculations with using quazipotential. But evidently potential was choosen 
unsuccessfully, as e turned out to be -0.9. Evidently this negative value of devoids it is 
clear physical sense, as a mixed parameters. Therefore this approach lost very much in 
heuristic understanding of obtained results. 

In Table 3 we present the results of hyperfine splitting calculation in heavy-quark 
systems. Namely exactly for these systems our Breit-Fermi approach must be true with 
maximum extent. The obtained results for the hyperfine splittings of the 5- wave states 
agree with measured splittings. As it is see most of our results have mainly predective 
character. In the Table 3 are shown that for 25'-states we obtain somewhat worse results 
for hyperfine splitting than in IS-states. This can be the result of a mixing of S and D 
waves. In 25* state, as it is shown in work [^], the mixing can give contribution of 10%, 
while in case of IS'-states the mixing correction is only of 1%. In other words we believe 
that taking into account the mixture of S and D waves, would considerable in proove our 
results. 

We suggest that the potential consist of a sum of vector and scalar parts. This idea 



of scalar- vector mixing was discussed in [0, 1T3-[|3]. The authors of these papers also 



came to the conclussion that its mixing parameter must be different from zero. Franzinis 



T7| show, that VCONF must be tottaly scalar, while VOGE must be tottaly vector, but 



nevertheless they cannot give an adeqyate description of data, concerning the fine split- 



ting. In series of papers |[T^-[^ Deoghuria and Chakrabarty had choosen the confining 
potential in the form 

We suggest that the potential consist of a sum of vector and scalar parts. This idea 



of scalar- vector mixing was discussed in [0, |T^-[^. The authors of these papers also 



came to the conclussion that its mixing parameter must be different from zero. Franzinis 
1% show, that VcoNF must be tottaly scalar, while Vqge must be tottaly vector, but 



nevertheless they cannot give an adeqyate description of data, concerning the fine split- 



ting. In series of papers |T^-[^ Deoghuria and Chakrabarty had choosen the confining 
potential in the form 



V = eVoGE+{l-e)Vc 



CONF 



and found e = 0.2, treating as an adjustable parameter. In this paper we have used the 
same approach for CJP-type potential and found e = 0.5. With this value of we described 
the hyperfine splitting of all mesons from heave to light ones. 
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